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Аннотация 

Введение. Существенная нелинейность химических реакций кислот и оснований определяет алгоритмы 
управления в режиме подкисления или подщелачивания, то есть периодической дозации минимального объема 
кислоты или щелочи. Такое регулирование может быть малоэффективным, т. е. допускает недостаточную или 
избыточную концентрацию контролируемого вещества. Статья посвящена проблеме точного регулирования 
водородного показателя в мини-биореакторах. Предлагается использовать цифровую модель системы 
управления кислотностью для подбора концентраций доливаемых растворов, определения методики 
регулирования и повышения точности. Цель работы — обеспечение требуемой точности регулирования рН 
в ш уйто мини-модели желудочно-кишечного тракта статического типа. 

Материалы и методы. Исходная структурная схема модели включает накопители и потоки. Она представляет 
собой базу для основных дифференциальных уравнений, характеризующих изменение объема и кислотности. 
Для корректировки показаний кислотности результирующего раствора по температуре создана статическая 
модель, основанная на полиномиальной аппроксимации экспериментальных данных методом наименьших 
квадратов. В прикладном пакете МаЧаб ЗипиПпК исследованы структурные элементы математической модели. 
На реальной системе ш уйго моделирования искусственного желудочно-кишечного тракта домашней птицы 
определены переходные характеристики для подтверждения адекватности математической модели. 
Результаты исследования. В рамках данной работы авторы создали и проанализировали нелинейную 
математическую модель изменения рН в биореакторе с учетом внешних управляющих воздействий. Потоки 
раствора соляной кислоты, раствора щелочи и слива раствора из реактора представлены как элементы 
дифференциального уравнения, описывающего накопление жидкости в реакторе. Для повышения точности 
решение доработали с учетом температурной зависимости водородного показателя. Предложена 
математическая модель дозирования на основе регулятора с каналами щелочи и кислоты. Полученные данные 
позволили генерировать объединенную модель процесса регулирования рН в биореакторе. Адекватность 
решения подтвердили опытным путем. В виде структурных схем показаны модели: регулятора РН, 
регулирования объема содержимого в реакторе и химических реакций. Сравниваются переходные процессы 
математической модели и реальной системы управления. Установлено, что переходные характеристики 
математической модели и реальной системы идентичны по времени регулирования. Относительная 
погрешность регулирования реальной системы составила 0,35 %, а математической модели — 0,1%, что 


соответствует требуемой точности регулирования -0,1 рН. В виде графиков показано влияние исследованных 
потоков на реакцию нейтрализации. 

Обсуждение и заключения. Предлагаемая математическая модель позволит подобрать оптимальные методы и 
алгоритмы регулирования кислотности, что ускорит создание регулятора нелинейного процесса регулирования 
водородного показателя. В будущем эти наработки можно интегрировать в комплексную цифровую модель 
всего искусственного желудочно-кишечного тракта домашней птицы для оптимизации алгоритмов управления 
(дозации, перемешивания, периодичности и т. п.), а также аппроксимации к объектам ш У1у0. 
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Введение. С начала 1980-х годов исследователи физиологии и функциональных особенностей внутренних 
органов человека пришли к пониманию, что фокус экспериментов будет смещаться от практик ш У1\о к Ш 
уго [1]. Актуальность разработки систем аппаратной имитации физиологических процессов обусловлена: 

— снижением затрат и времени на проведение исследований: 

— возможностью единовременного проведения однотипных экспериментов; 

— акцентом на этичность исследований, что предполагает ограничение экспериментов на живых существах. 

К началу ХХПГвека сформировалось научное направление по разработке методик и оборудования для 
моделирования и исследований т уго [1]. 

В настоящее время известны статические и динамические ш уго модели желудочно-кишечного 
тракта (ЖКТ). В первом случае речь идет о моделировании процессов пищеварения при использовании жидких 
питательных сред [2]. Динамические модели, в отличие от статических, позволяют моделировать 
перистальтические движения в ЖКТ и используются при работе с пищевыми матрицами в виде суспензий, 
включающих твердые частицы и волокна. 

Цель данной работы — обеспечение требуемой точности регулирования рН в ш \о мини-модели ЖКТ 
статического типа. 

Для достижения цели решаются следующие задачи: 

— разработка нелинейной математической модели изменения рН в биореакторе с учетом внешних 
управляющих воздействий; 

— исследование математической модели методами имитационного моделирования; 

— доработка математической модели для повышения ее точности путем учета температурной зависимости 
водородного показателя; 

— создание математической модели процесса дозирования на основе двухканального конкурентного 
регулятора (по каналам щелочи и кислоты); 

— разработка объединенной математической модели процесса регулирования рН в биореакторе; 

— проведение экспериментов для подтверждения адекватности объединенной математической модели и 
определения достигнутой точности регулирования. 

Материалы и методы. С точки зрения моделируемого объема модели ЖКТ подразделяются на микро- (до 
50 мл), мини- (от50 до 400 мл) и макросистемы (свыше 400 мл) [3]. Отдельно следует упомянуть 
микрофлюидные системы [4], но они выходят за рамки данной работы. 

Патентный анализ позволяет назвать первое зарегистрированное изобретение в данной сфере — «т уйго 
модель ш У\о пищеварительного тракта»'. Эта динамическая система предшествовала появлению систем типа 
ТИМ? [5]. Идея развивалась в моделях ООМ [6], ТИМ? [5], НО$ [7] и др. Статические модели (ОТООТ [8], 
У1МОТ [9], 5НИМЕ [10], АВСОЕ. [11] и др.) ушли от сложного моделирования перистальтики и сфокусировались 
на более точном моделировании собственно процессов пищеварения. Первая причина — техническая 
сложность и дороговизна реализации систем с динамическим сжатием содержимого. Вторая — развитие 
методов математического моделирования и микроконтроллерных систем управления, позволивших мониторить 
среду в пробирке (реакторе) с высокой точностью и стабильностью [4]. 

В данной работе рассматривается решение частной задачи математического моделирования процесса 
регулирования рН для ш УЙго мини-модели желудочно-кишечного тракта (ЖКТ) статического типа [12]. 
Точность управления кислотностью химуса определяется жесткими требованиями к условиям экспериментов п 
уго. Точность автоматического регулирования рН в биореакторе не может быть ниже +0,1 единицы.” Это 
принципиально для правдоподобия условий процессов пищеварения и соблюдения требований к 
повторяемости результатов опытов. При этом рН может быть от слабощелочной (7,9 рН, ротовая полость) до 
сильно кислотной (1,3 рН, желудок). 


Для моделирования реакции нейтрализации воспользуемся методом накопителей и потоков [13]. Его база — 
это законы сохранения (накопители), а также математическое описание закономерностей, изменяющих 
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количество накапливаемого признака (потоки). Структурная схема данной модели с двумя накопителями 
представлена на рис. 1. 


(1) Объем 
кислоты 5 
1 
Емкость (3) Сброс 
(3) Объем 
щелочи 
Реактор 
(4) Поток 
Маон о, (6) Процессы 
Раствор 
(5) Поток пищеварения 
НС1 


Рис. 1. Структурная схема модели с двумя накопителями 


Реактор как физический объект накапливает некоторое количество: 
— жидкости (01); 
— ионов водорода (О). 
В этой емкости идет контролируемый процесс пищеварения с изменением рН и объема жидкости. 
На первый накопитель действуют три потока: 
— 1-й — объем доливаемой кислоты; 
— 2-й — объем доливаемой щелочи; 
— 3-й — объем, выливаемый из реактора. 
На второй накопитель также влияют три потока: 
— 4-й — подача насосом гидроксида натрия [2]; 
— 5-й — поступление соляной кислоты [2]; 
— 6-й — изменение рН в процессе пищеварения. 
Последний поток представляет собой возмущающее воздействие и не исследуется в данной работе. 
Результаты исследования. Рассмотрим математическую реализацию модели для описания физико- 
химических процессов. 
Объем жидкости в реакторе: 
У=5. 1. (1) 
Изменение объема жидкости в реакторе: 
аУ =5 -ай. (2) 
Для удобства назовем: 
— поток раствора соляной кислоты — 41 = С1; 
— поток раствора щелочи (гидроксид натрия) — 42 = Со; 
— поток слива раствора из реактора — 42 = Сз. 
Таким образом, дифференциальное уравнение, описывающее накопление жидкости в реакторе: 
оао (3) 
% 5 
Для определения динамики изменения концентрации раствора в реакторе необходимо найти дифференциал 
количества вещества: 


УЕУ.С, (4) 
где У — это объем вещества, л; С — молярная концентрация, моль/л. 
Из (4) и (1) получаем: 
аа ь: | (5) 
& Ш ф @ 


На изменение концентрации раствора влияют потоки щелочи и кислоты с определенной концентрацией. 
Значит: 


ай ас 
5-м. | - (@-Сыи-,-Сл). (6) 
Из (6) выведем изменение концентрации и упростим выражение: 
ас 1 
=— (С, (Сион ие (п 
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Представим систему дифференциальных уравнений, описывающих динамику изменения объема и 
концентрации раствора в биореакторе: 


Ца +6,-С,), 
@ 5 
(8) 
ас _ 1 
и? - яр (С (мон -С+ „(С — Сна )). 


Зная основы химии и принимая во внимание (8), можно определить концентрацию катионов Н ([Н*]). 
Концентрации ионов водорода ([Н*]) и гидроксид-ионов ([ОН ]) в дистиллированной воде при 25 °С равны и 
составляют по 10” молыл. Это следует из определения ионного произведения воды [11]: 


к = [Н*|х [ОН ] =Ю “тор /Р (а 25 °С). (9) 
Также стоит учитывать, что: 
Т=[ОН`]-[Н"] = [Ма] СГ]. (10) 
Теперь найдем [Н"] (обозначим, как Си) через квадратный корень из предыдущих выражений с учетом 
концентрации раствора по уравнению (4) [14]: 


О и 
С, =УК,Т=0, (11) 
с, =-5 т 


Здесь ионы 1= 4. 

Для создания системы регулирования с обратной связью по кислотности в виде водородного 
показателя (рН) используем отрицательный десятичный логарифм из концентрации, рассчитанной по 
уравнениям (7), и получим кислотность раствора при нормальных условиях, то есть при 25 °С. 

Однако содержимое биореактора не всегда находится в нормальных условиях. Например, для создания 
искусственной среды, свойственной ш У1уо среде желудочно-кишечного тракта домашней птицы, в реакторе 
поддерживается температура 42 °С. Поэтому на высокоточных измерителях рН получили матрицу изменения 
рН буферных растворов в зависимости от температуры. Показания фиксировали по измерительной системе с 
термокомпенсацией для применяемого электрода [13]. С этой целью использовали данные эксперимента и 
сведения из открытых источников°. 

В итоге построили статическую модель. Ее получили полиномиальной аппроксимацией методом 
наименьших квадратов. Зависимость рН от температуры визуализирована на рис. 2 [14—16]. 


Реальный рН > с 


Рис. 2. Поверхность, описывающая зависимость рН от температуры (12) 


3 Ашотанс Тетрегацие Сотрепзаноп ш РН Меазигетеи/ Нопфа// Вопфа.сот: [сайт]. ВТ: Вирз:/Амум.Вофа.сот/езр/мацег- 


диаНеу/зирронАесьщса]-Нрз/БепсН-пеегз/ащюотайс-{етрегаиге-сотрепзаноп-ш-рЬ- 
шеазигетеп И? ии _зоигсе=ив\у&ини_тедат=301&хани_сатраеи=иб\-ге@ тес! (дата обращения: 25.12.2022). 
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Е(®) = р00+ р10-х+ р01- у+ р20-х? + р. х. у+ р02- у? + 2) 
+р30-х* + р21-х’-у+р12-х- у”. 


Коэффициенты уравнения представлены в таблице 1. 


Таблица 1 
Коэффициенты уравнения (11) 
Коэффициент Значение 
200 0,154 
р10 —0,00511 
ро1 0,899 
р20 —2,998 10° 
РИ 0,00333 
р02 0,0130 
р30 -3,226. 10” 
р21 1,588 . 10° 
р12 —4,8. 10“ 


Для определения переходных процессов в системе построим математическую модель в прикладном пакете 
МаЧаь ЗитиНак (рис. 3): 


ОН 
ОН 
Реальный РН У РН(25С) 
Н 
Модель системы Модель Н 
управления дозацией изменения 


объема 


РН Реальный рН 


Статическая модель 
зависимости рН от 
температуры 


Модель 
химических реакций 


Рис. 3. Общая структура математической модели в МаНаб ЗипаПаК 


Рассмотрим созданную систему по блокам. Первый — блок двухканального регулятора по кислотности с 
обратной связью. В цепь обратной связи введены шумы для моделирования случайных событий реальной 


системы измерения кислотности и погрешности измерительных приборов. 


ила сы 
Требуемый рН ПП 
нь 


Контроллер Н 


+ ий Ошибка считывания рН 
Реальный рН 
ь Со 


Рис. 4. Структура регулятора рН, блок «Модель системы управления дозацией» 


Второй блок — это система контроля объема раствора в реакторе. Функция Уоате свапее то4е 
моделирует включение откачивания. Так обеспечивается постепенный сброс раствора до необходимого объема 


при его приближении к максимуму. 
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Начальный 
объем 


Контроллер сброса 


Уо 


у ощ 


х 1/5 (реактор) 


Объем в высоту 
столба жидкости 


Графики 


и 


У 


Рис. 5. Структура модели регулирования объема содержимого в реакторе, блок «Модель изменения объема» 


Третий блок — основной. Он реализует расчет химической реакции кислоты и щелочи при добавлении в 
реактор и вычисление рН конечного раствора при нормальных условиях. Функция МеибаПтайоп определяет 
концентрацию ионов водорода в растворе [12]. Блок Сопсепиайоп ш рН преобразует концентрацию в 
водородный показатель рН. 


-(Расчет реакции) 


Концентрация рН 
при н.у. 


Нео 


РН (25 С) 


Рис. 6. Структура модели, описывающей химические реакции, блок «Модель химических реакций» 


Функция статистической модели в четвертом блоке (Тетрегаеиге соггесНоп) позволяет учесть температурные 


изменения раствора. 


Температура 
раствора 


Температурная коррекция 
водородного показателя 


Реальный РН 
раствора 


Рис. 7. Структура модели, описывающей изменение кислотности от температуры, блок 


«Статическая модель зависимости РН от температуры» 
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Для определения адекватности переходных характеристик результирующей математической модели 
проводились эксперименты на имитационной модели. Реальная модель представлена системой 
автоматизированного управления искусственным ЖКТ домашней птицы с одним блоком управления и 
биореактором. Оборудование представлено на рис. 8. 


а) 6) 


Рис. 8. Оборудование для проведения исследований: а — система управления параметрами среды в биореакторе; 


б — элементы биореактора. 1— управляемые перемешивающие устройства; 2— нагреватели; 
3— порты для датчиков и системы дозации жидкостей 


Опишем методику проведения эксперимента. В математической модели учитываются характеристики 
реальной системы: 

— заданная кислотность; 

— параметры реактора (площадь основания и высота); 

— концентрации доливаемых веществ (кислоты и щелочи); 

— начальная концентрация раствора; 

— начальный объем раствора; 

— температура раствора на протяжении всего эксперимента; 

— дискретность процесса регулирования [14]. 

Математическая модель представляет собой идеальную систему: смешивание растворов в реакторе 
происходит мгновенно и равномерно. В реальных условиях это невозможно, поэтому для равномерного 
растворения добавленной кислоты или щелочи необходимо принудительное перемешивание. В имитационной 
модели данная проблема решается магнитной мешалкой с управлением $САПА“. 

К сожалению, такой тип перемешивающего устройства наводит на рН-электрод электромагнитные помехи. 
Для исключения помех устройство перемешивания кратковременно отключалось. Дискретность измерений и 
работы регулятора составила 20 секунд, что обусловлено экспериментально установленным временем 
нормализации показаний рН электрода в реакторе (по документации электрода возможна дискретизация до 1- 
2 минут) ”°. В целях обеспечения безопасности регулятор реальной системы ограничен в максимальном объеме 
дозируемой кислоты или щелочи, что учитывалось в модели. 

На рис. 9 сравниваются переходный процесс математической модели и реальной системы управления. РН- 
электрод — это чувствительный элемент, на который воздействуют внешние электромагнитные помехи. 


4 От англ. Зирегу1зогу Сопто! Апа Рава Асаш$#юоп — диспетчерское управление и сбор данных. Программный пакет для создания или 
поддержки работы систем сбора, обработки, отображения и хранения информации об объекте мониторинга или управления. 

° РЬ-4502< РН шеег саНьаноп поез / Шопе01,Ы0е / И опе01510°. Му @есиотсз Фу ргодесз: [сайт]. ОВ: 
Брз://Фоп201 .6102/2019/04/26/рь-4502с-рЬ-теег-саШгаНоп-побе5/ (дата обращения: 30.11.2022). 

° Н7Отм рН Сисай ПаазВее / АЧаз Заепийс МС // НШез.аМаз-зслепайс.сот: [сайт]. ОВГ:  Бирз://ез.аЧаз- 
зс1епайс.сош/рН_ЕЙО_РаазВееЕ.рАЕ (дата обращения: 30.11.2022). 
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Рис. 9. Переходные характеристики математической модели и реальной системы при регулировании от 6,2 до 2,15 рН 


Один из экспериментов проводился при заданной кислотности раствора 2,15+0,05 рН. Кислотность 
исходного раствора — 6,2 рН. Температура раствора в реакторе — 22,25 °С. Для изменения кислотности 


использовались растворы 


0,1 моль/литр МаОН и 0,1 моль/литр НС. 


дискретном пропорциональном регулировании. 


Алгоритм дозации заключался в 


На рис. 10 представлены интегралы по потокам из дифференциального уравнения (7). 


6:10-8 


4-10 


Интегралы по потокам 
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300 


350 


400 


450 — Время, с 


Рис. 10. Влияние каждого потока на реакцию нейтрализации в ходе экспериментов 


Обсуждение и заключения. Анализ результатов моделирования показал, что потоки с концентрацией 


раствора в реакторе не вносят весомого вклада в регулирование всей системы, поэтому ими можно пренебречь. 
Тогда система дифференциальных уравнений цифровой модели (8) приобретет следующий вид: 


а<_ 1 


а 1 
в = < (6, +С(, —6. }, 
(С, } Стон —С, "Сна ). 


@& 51 


(13) 


По результатам экспериментов с учетом концентрации раствора в реакторе относительная погрешность 


установившейся кислотности раствора для математической модели составила 0,1 %, а для реальной системы с 
тем же алгоритмом — 0,35 %. 


Итак, в математической модели учтем погрешность показаний электрода и пренебрежем перемешиванием 


растворов. Абсолютная погрешность результатов регулирования находится в Допустимых пределах и не 


превышает 0,1 РН от заданной величины. Это означает, что разработанная математическая модель адекватно 


описывает кривую реакции нейтрализации и в будущем позволит подобрать оптимальные алгоритмы 
регулирования, что было бы трудоемким процессом в реальной системе. 
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Разработанная модель решает поставленную задачу. Она представляет собой математический инструмент 
для создания методик и алгоритмов регулирования РН с учетом температуры содержимого мини-реактора, 
начальной концентрации основного раствора и концентраций доливаемых растворов. Такой подход позволяет 
совершенствовать существующие алгоритмы управления. 
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